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Abstract. Using XMM-Newton observations, we investigate the scaling and structural properties of the ICM entropy in a sample 
of 10 nearby (z < 0.2) relaxed galaxy clusters in the temperature range 2-9 keV. We derive the local entropy-temperature 
(S—T) relation at R = 0.1,0.2,0.3 and 0.5/?200- The logarithmic slope of the relation is the same within the lcr error at all 
scaled radii. However, the intrinsic dispersion about the best fitting relation is significantly higher at 0. IRioo- The slope is 
0.64 ± 0. 1 1 at 0.3 R200, in excellent agreement with previous work. We also investigate the entropy-mass relation at density 
contrasts d - 5000, 2500 and 1000. We find a shallower slope than that expected in simple self-similar models, which is 
in agreement with the observed empirically-determined entropy-temperature and mass-temperature scaling. The dispersion 
is smaller than for the S-T relation. Once scaled appropriately, the entropy profiles appear similar beyond ~ O.I.R200, with an 
intrinsic dispersion of ~ 15 per cent and a shape consistent with gravitational heating (S (r) 2 r 1 1 ). However, the scatter in scaled 
entropy profiles increases with smaller scaled radius, to more than 60 per cent at R < 0.05i?200- Our results are in qualitative 
agreement with models which boost entropy production at the accretion shock. However, localised entropy modification may 
be needed to explain the dispersion in the inner regions. 
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1. Introduction 

X-ray observations of the hot, gaseous intracluster medium 
(ICM) have been telling us for well over a decade that physical 
processes other than gravity are acting to modify the properties 
of the cluster population. An understanding of the source(s) of 
this modification is of great importance for our understanding 
of cluster formation and evolution, and is essential for the use 
of clusters as precision cosmological probes. 

Entropy is important because, together with the shape of the 
potential well, it is the quantity which dictate s the observ ed X- 
ray properties of the ICM in galaxy clusters dVoitl I2005). The 
intracluster entropy is generated in shocks a s gas is drawn into 
the gravitational potential of the cluster halo (jTozzi & Norman, 
200 UlBorgani et al l bOOU IVoit et all I20021 l2003U VoilU2005l 
Bor gani et all l2005h . thus it reflects the accretion history of 
the ICM. However, the entropy distribution also preserves key 
information regarding the influence of non-gravitational pro- 
cesses. 

X-ray measurements of the entropy 1 at O.I-R200 (hereafter 
S 0.1) showed that the entropy of the coolest systems is con- 
siderably higher than that available from gravitational col- 
Send offprint requests to: G. W. Pratt, e-mail: gwp@mpe.mpg . de 

1 Throughout this paper the entropy is defined as 5 = kTn^ 2/3 . This 
quantity is related to the true thermodynamic entropy by a logarithm 
and an additive constant. 



lapse JPonman et allll999tlLlovd-Davies etaiT.120001) . and that 
the entr opy-temperature (S- T) relation is shallower than ex- 
pected dPonman et all [2003). More recent spatially resolved 
entropy profiles indicate that the entropy is higher throughout 
the ICM, and that, outs ide the core regions, entropy profiles 
are structurally similar JPo nman et all l2003t IVoit & Ponmanl 
l2003HPratt & ArnaudlEo03U2005HPiffareti et alll2005tl . At the 
same time the scatt er in So 1 at a given temperature can be up 
to a factor of three JPonman et alll2003l) . 

Possible entropy modification mechanisms have histori- 
cally included preheating, where the gas has been heated be- 
fore being accreted into the p otential well, b y early supernovae 
and/or AGN acti vity (e.g.. iKaisert 1199 ll |Eyrard & Henrvi 
ll99lUValageas & Silklil999l) . internal heating a fter accretion 
(e.g.JMetzler & Evrardll 19941) . and cooling (e.g. JPearce et all 
2000). The lack of isentropic core entropy profiles in groups 
and poor clusters has shown that simple preheat ing is unli kely 
to be the sole explanation of the observations (Pon man et all 
l2003UPratt & Arnaudll2003Ll2005l) . Since cooling-only models 
gene rally predict a higher st ellar mass fraction than observed 
(e.g.. lMuanwong et all l2002t) . attention is now focussing on the 
interplay between cooling and feedback. Further high quality 
observations are needed in order to distinguish between these 
different entropy modification mechanisms. 

The present paper is the third in a series based on 
XMM-Newton observations of ten clusters in the temperature 
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Fig. 1. Cluster entropy profiles obtained from the deprojected, PSF corrected temperature profiles and the best fitting analytical model for the 
gas density. Solid lines, included to improve visibility, are entropy profiles obtained from analytic model fits to the temperature and density 
information. On the left, the profiles are plotted in physical units QCj^ kpc); on the right, they are plotted in units of ^200 • 



Table 1. Basic cluster data. 



Cluster 


z 


T 

(keV) 




A1983 


0.0442 


2.18 ± 


0.09 


A2717 


0.0498 


2.56 ± 


0.06 


MKW9 


0.0382 


2.43 ± 


0.24 


A1991 


0.0586 


2.71 ± 


0.07 


A2597 


0.0852 


3.67 ± 


0.09 


A 1068 


0.1375 


4.67 ± 


0.11 


A1413 


0.1427 


6.62 ± 


0.14 


A478 


0.0881 


7.05 ±0.12 


PKS0745 


0.1028 


7.97 ± 


0.28 


A2204 


0.1523 


8.26 ± 


0.22 



range from 2 keV to 9 keV. In Pointe couteau. Arnaud & Pratt 
(2005, hereafter Paper I), we investigated the shape and 
scaling properties of the dark matt er distribution, and in 
lArnaud. Pointecouteau & Prattl d2005l hereafter Paper II) we 
examined the relation between the total mass and the X- 
ray te mperature. In th e present paper we extend the work of 
|Pratt& Arnau l ll2004 by exploring the entropy scaling prop- 
erties over a larger temperature/mass range. 

All results given below were calculated assuming a ACDM 
cosmology with Q m = 0.3 and Q.\ = 0.7 and Hq = 
70 km s _1 Mpc~'. Unless otherwise stated, errors are given at 
the 68 per cent confidence level. 

2. Sample and analysis 

The sample (Tabled comprises 10 systems ranging in tem- 
perature from 2 keV to ~ 8.5 keV. In Paper I, spatially re- 
solved temperature and density data were used to calculate 



gravitating mass profiles. These were which then fitted with 
an NFW model, yielding values for M200 and R200 ('virial' 
mass and radius corresponding to a density contrast 2 of 200). 
The mean temperature of each system, r spe c, needed to in- 
vestigate the scaling properties of the profiles, was estimated 
from fits to a spectrum built from all events in the region 
0.1/?200 < r < 0.5/? 200 (see the discussion in Paper II). The 
data reduction steps are described in detail in Paper I. 

We calculated the entropy profile of each system using the 
deprojected, PSF corrected temperature profile and an analyti- 
cal model for the gas density derived in Paper I. Analytic gas 
density model parameters are listed in Appendix A. In Fig ^ 
the raw entropy profiles are shown in physical units (hri kpc), 
and in terms of the measured virial radius, ^200- All profiles 
increase monotonically with radius and, while the slope of the 
profile becomes shallower towards the centre in some of the 
clusters, none shows an isentropic core. 

3. Entropy scaling properties 

3.1. Introduction 

In a self-similar cluster population expected from simple grav- 
itational collapse, the scaled profiles of any physical quantity 
coincide, and thus measures of these quantities at any scaled 
radius should correlate with global quantities such as the mean 
temperature, or the virial mass. In this standard model, the en- 
tropy is expected to scale as S oc h{z)~~ A ^T or S oc h{z)~ 2 ^M 2 ^ . 

In purely theoretical terms, the most fundamental char- 
acteristic of a cluster is its mass. Previous investigations 

2 The density contrast, 6 = 3M(< r)/4nr 3 p c (z), where p c (z) = 
h 2 (z)3H 2 /S7TG, and h 2 (z) = n M (l + z) 3 + fi A . 
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Table 2. The S-T relation. Data were fitted with a power-law of the 
form h(z) 413 S = A x (kr/5keV)", where kT is the overall spectro- 
scopic temperature in the 0. 1 - 0.5 R200 region. Errors in entropy and 
temperature are taken into account. Results are given for the WLS, 
BCES, and WLSS regression methods (see text). Raw, statistical and 
intrinsic scatter about the best fitting relation in the log-log plane are 
given in the last columns; the numbers in parentheses indicate the per- 
centage intrinsic scatter. 



Radius 


A 






a 




0" log 




-^200 


keV cm 2 






raw 


stat 


int 


WLS 
















0.1 


227 


±5 


0.58 


±0.05 


0.079 


0.030 


0.073 


0.2 


478 ± 


12 


0.73 


±0.06 


0.058 


0.035 


0.047 


0.3 


770 ± 


24 


0.71 


±0.07 


0.074 


0.043 


0.060 


0.5 


1510 ± 


90 


0.68 


±0.12 


0.070 


0.078 




BCES 
















0.1 


230 ± 


17 


0.49 


±0.15 


0.082 


0.030 


0.076 (19%) 


0.2 


485 ± 


22 


0.62 


±0.11 


0.063 


0.034 


0.052 (13%) 


0.3 


798 ± 


44 


0.64 


±0.11 


0.078 


0.043 


0.065 (i6%) 


0.5 


1560 ± 


83 


0.62 


±0.08 


0.074 


0.078 




WLSS 
















0.1 


229 ± 


16 


0.47 


±0.14 


0.083 


0.030 


0.077 


0.2 


480 ± 


23 


0.67 


±0.10 


0.059 


0.035 


0.048 


0.3 


786 ± 


43 


0.69 


±0.12 


0.075 


0.043 


0.061 


0.5 


1510 ± 


90 


0.68 


±0.12 


0.070 


0.078 





of the entropy scaling relations were hampered in this re- 
spect because they did not have accurate mass information. 
The present observations, which have excellent mass data 
(Papers I and II), thus represent an ideal opportunity to ex- 
amine the entropy-mass relation. In addition, the spatial res- 
olution and radial coverage in the entropy profiles are signif- 
icantly i mproved, as c ompared with previous ROSAT/ASCA 
studies ( Pon man et all Eo03l) . In this Section, we investigate 
both the entropy-temperature (h(z) 4 ^S-T spec ) and the entropy- 
mass (h(z) 2 ^S-M2oo) relations. We derive the slope and nor- 
malisation of each relation, at various scaled radius/density 
contrasts. We also estimate the raw, statistical and intrinsic 
scatter about the relations. 

3.2. Method 

The entropy estimates at a given scaled radius (x = r/7?2oo) or 
density contrast (6) have to be calculated taking into account 
the uncertainties on the temperature, the analytical gas density 
model, and the measurement of the scaled radius/density con- 
trast itself. We thus estimated the entropy at a given x or 6, 
and the associated errors, using a Monte Carlo method. The 
scaled radius r(x)/Rg values and corresponding uncertainties 
were estimated from NFW-type modelling of the mass profile 
(see Papers I and II for details), and subsequently randomised 
assuming a Gaussian distribution with sigma equal to the lcr 
error. The temperature profile was then similarly randomised 
taking into account the observed uncertainties, and the tem- 
perature estimated at the relevant (randomised) r(x)/R$ using 
spline interpolation. The entropy was then calculated in the 
usual manner, S — kT/n e , including an additional 5 per cent 
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Fig. 2. The S-T relation measured from a sample of 10 clusters cov- 
ering a temperature range from 2 to 9 keV. The S-T relation is shown 
for different fractions of R 2 oo- Measurements are plotted with error 
bars. At each radius, the best-fitting power-law relation, derived tak- 
ing account the errors in entropy and temperature, is overplotted. Solid 
line: BCES regression method; dashed line: WLS method; slopes and 
intercepts are given in Table|2| 

systematic error on n e (corresponding to the uncertainties in 
the surface brightness profile modelling; see Paper I). This pro- 
cedure was undertaken 1000 times for each cluster; the final 
entropy value at each xj 6 was calculated from the mean and 
standard deviation of the randomised values at that radius. 

The temperature profiles of all systems except A1983 and 
MKW9 are detected up to ~ 0.5 R200 (the scaled radius/density 
contrast of maximum detection of A1983 and MKW9 are 
0.38^2oo/1455(5 and 0.41/? 2 oo/ 1401(5, respectively). The tem- 
perature profiles of these clusters had thus to be extrapo- 
lated to Q.5R 200/1 0005. To avoid spuriously high/low temper- 
atures, we only allowed extrapolations within +1.5 times the 
mean temperature of the cluster (well with in the range of tem- 
perature gradients ob served to date, e.g.. lArnaud et alj [2005: 
IVikhlinin et all2005t) . 

3.3. The S-T relation 

Our data have the required temperature coverage to investigate 
the entropy-temperature relation. For easier comparison with 
previous work, we first investigated the S-T relation at differ- 
ent fractions of the virial radius, viz., 0.1, 0.2, 0.3, and 0.5 ^200- 
We fitted a power law model of the form 

/z(z) 4/3 S x =A[r spec /5keVr (1) 

to the data, r spec being defined as in Sect. |2] We chose a 
pivot point of 5 keV for consistency with the M—T relation of 
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lArnaud et alJj2005h . The fit was performed using linear regres- 
sion in the log-log plane, taking into account the errors on both 
variables. We have also computed the raw and intrinsic scatter 
about the best fitting relations in the log-log plane. To estimate 
the raw scatter, we used the vertical distances to the regression 
line, weighted by the error. The intrinsic scatter was computed 
from the quadratic difference between the raw scatter and the 
scatter expected from the statistical errors. The resulting values 
are given in Table|2] 

In undertaking the fits, we first used the classical weighted 
least square method, WLS (y 2 estimator) impl emented in 
the th e routine FITEXY from numerical recipes dPress et all 
Il992l) . This regression method is strictly valid only if the in- 
trinsic scatter is negligible as compared to the statistical scat- 
ter. This is generally not the case (see Tabled . We thus also 
consi dered the orthogonal BCES method ( Akritas & Bershadv, 
Il996l) . with bootstrap resampling. This is the least-biased re- 
gression method when both measurement errors and intrinsic 
scatter are present. However, no error weighting is performed 
on individual data points. We thus also considered a variation 
of the WLS method, where a constant term, corresponding to 
the intrinsic scatter, is added quadratically to the log(5 ) error 
(WLSS method). The best fit is determined such that the re- 
duced^ 2 , after minimization over A and a, is^ 2 ed = 1. 

The best fitting slopes and intercepts for the different meth- 
ods are listed in Table|2] Figure|2]shows the data and the best- 
fitting power law relation for each radius under consideration, 
obtained using the standard WLS and BCES methods. The 
three methods give best fitting parameters consistent within 
the lcr errors. However, the errors are underestimated with the 
WLS method, especially at low radii. This is a consequence of 
the high x 2 value obtained, reflecting the high intrinsic scatter 
at small scaled radius (e.g. x 1 - 54 for 8 d.o.f at 0.1 /?20o)- In 
this case the standard A^ 2 = 1 criteria to derive parameter er- 
rors is not valid. The BCES and WLSS give results in excellent 
agreement with each other, both in terms of best-fitting values 
and associated uncertainties. In the following, for easier com- 
parison with previous work, we will refer to results obtained 
using the commonly used BCES method. 

The slope of the entropy-temperature relation is incompati- 
ble with the standard self-similar prediction at all radii at w hich 
we ha ve measured it, confirming the results of IPonman et alJ 
J2003I) . Our best-fitting S QA -T slope (a _ 0.49 ± 0.15) is 
slightly shallower than that found by IPonman et alJ from un- 
weighted orthogonal regression on individual data points (a = 
0.65 + 0.05). A better agreement is observed with the slope de- 
rived by these authors from data grouped in temperature bins 
(a = 0.57 + 0.04). However, the slopes are consistent once un- 
certainties are taken into account. 

Beyond r > 0.2 R 200, the slope remains remarkably stable. 
At 0.3/?200, our gas density and temperature measurements are 
well constrained, and do not suffer potential systematic prob- 
lems connected to the correction for PSF and projection effects. 
At this radius, we are well outside the strong cooling cores 
found in the majority of the objects in our sample. In addition, 
we note that there are no significant temperature gradients at 
thi s radius in the h i gh re solution temperature profiles produced 
bv lVikhlinin et al.l J2005ll from Chandra observations. For this 



Table 3. The S-M relation. Data were fitted with a power-law model 
of the form h(z) 2n S 6 = B 6 x (M 200 / 5.3 x 10 14 M Q ^, where M 200 is 
the total mass obtained from NFW model fits to the measured mass 
profiles. Best-fitting relations are given for the BCES estimator, taking 
into account errors in mass and entropy. Results are given for the full 
sample, and for the sample excluding MKW9, which is a significant 
outlier. The numbers in parenthesis indicate the percentage intrinsic 
scatter. 



6 


B 






O" log 




keV cm -2 




raw 


stat int 



Full 

sample 

"3000 471 ± 18 0.36 ±0.10 0.058 0.034 0.046 <n%) 

2500 765 ±30 0.37 ±0.10 0.059 0.041 0.043(10%) 

1000 1460 ±47 0.36 ±0.06 0.059 0.065 - 
Excl. 
MKW9 

~5000 459 ±37 0.43 ±0.10 0.039 0.035 0.017(4%) 

2500 741 ± 43 0.44 ± 0.08 0.035 0.041 - 

1000 1430 ± 46 0.41 ± 0.04 0.042 0.063 - 



reason we consider the measurements at 0.3/? 200 to be the most 
reliable. The slope of the S0.3-T relation, obtained using the 
BCES method , is Sm oc x 0M±0Al , in excellent agreement with 
that found by Ponma n et alJ J2003I) . 

Figure [5] shows that there is noticeable scatter about the 
S0.1—T relation. Table[2]makes clear that the scatter is reduced 
at larger scaled radius. The intrinsic scatter remains the dom- 
inant contributor to the dispersion in all relations, except at 
0.5 /?20o- 

3.4. The S -M relation 

In investigating the S-M relation, it is more logical to work 
in terms of density contrast 5. We used 6 = 5000, 2500, and 
1000, which correspond to average fractions of 0.20+0.01 Rim, 
0.29 + 0.02 #200 and 0.47 ± 0.02 fl 2 oo for the present sample, 
and fitted a power law of the form 



h(z) Vi S s = B s [M 200 /53 x 1O 14 M } 



(2) 



to the data. The pivot point of 5.3 x 10 14 M© corresponds to 
a temperature of 5 k eV using the best-fitting M-T relation of 
lArnaud et alJ i2005h . In FigE] we show the data and the best- 
fitting power law relation for each density contrast under con- 
sideration. Since MKW9 appears to be a clear outlier, we also 
fitted the relation excluding this cluster. The best fitting slopes 
and intercepts for the S-M relation, with and without MKW9, 
are listed in Tabled 

The best- fitting S-M relation is shallower than the predic- 
tion from standard self-similar models. The results for the full 
sample are entirely consistent with the observed S -T and M-T 
relations (Sect. [3~31 and Paper II). Excluding MKW9, the slope 
is slightly, although not significantly, steeper. 

As can clearly be seen in Fig. [3] and as quantified in Table|3j 
the intrinsic dispersion about the S -M relations is less than that 
for the S-T relations. This is particularly true for the sample 
excluding MKW9. This is not an artifact of the larger statistical 
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Fig. 3. The S-M relation measured at various density contrasts. M200 
values were measured from NFW model fits to the measured mass pro- 
files (Papers I and II). Symbols with error bars represent the measure- 
ments. Lines show the best-fitting power-law S-M relation: solid: full 
sample; dashed: excluding MKW9. Slopes and intercepts are given in 
Tabled 
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Fig. 4. Scaled entropy profiles. The radius is scaled to i?200 measured 
from the best-fitting mass models (Papers I and II). The entropy is 
scaled using the empirical entropy scaling S oc h(zy 4/3 T°g 5 , using 
the global temperature, T, as listed in Table Q in units of 10 keV. 
The shaded grey area corresponds to the region enclosed by the mean 
plus/minus the lcr standard deviation. The dashed line denotes 5 oc 
R im . 



mass errors compared to those on the temperature: the disper- 
sion is similarly smaller for the raw scatter. Table [5] also indi- 
cates that, excluding MKW9, the dispersion is now dominated 
by the statistical scatter at all density contrasts. 

4. Scaled entropy profiles 

We have confirmed that the S-T relation is shallower than ex- 
pected in the standard self-similar model, and we have shown 
that this is also true for the S-M relation. In this Section, we 
examine the structural properties of the entropy profiles. We do 
this by looking at the entropy profiles once scaled by the best 
fitting relations. 

4.1. Scaled entropy versus r//? 2 oo 

In Fig.|4]we show the profiles scaled using the relation S oc 
h{zY 4 ^ T®q 5 , where T\q is the global temperature measured 
in units of 10 keV. This relation is consistent with our data 
(Table and allows us to compare our results with previ- 
ous work. As an initial measure of the scatter in scaled en- 
tropy profiles, we estimated the dispersion at various radii in 
the range 0.02-0.45 /?200- The shaded area in Fig @] shows the 
region enclosed by the mean plus/minus the lcr standard devia- 
tion. Figure@]shows that, outside the core regions, the entropy 
profiles present a high degree of self-similarity. The relative 
dispersion in scaled profiles remains approximately constant at 
< 20 per cent for r > 0. 1 ^200, in excellent agreement with 
the dispersion found in a smaller subsample bv lPratt & Ar naud 



d2005l) . In the core regions, however, the dispersion increases 
with decreasing radius to reach > 60 per cent at ~ 0.02 /?20o- 

We next fitted the scaled profiles with a power law using 
the BCES method. The fit was performed in the log-log plane 
taking into account the errors on both variables. Fitting the data 
in the radial range r > 0.01 /?200 we find a slope of 1.08 + 
0.04, with a dispersion of ~ 30 per cent about the best fitting 
line. The slope is not significantly changed (1.14 + 0.06) if the 
data are fitted in the radial range r > 0.1 R200, but the intrinsic 
dispersion is two times smaller (14 per cent). 

For comparison, IPratt & Arnaudl J2005h . using a subsam- 
ple of the current data (A1983, MKW9, A271 7, A1991 and 
A1413 ), found a slope of 0.94 + 0.14, while IPiffareti et all 
d2005l) . using the same scaling, find a slope of 0.95 + 0.02 
for their sample of 13 cool core cluste rs. In contrast, broken 
power-law behaviour has been found bv lMahdavi etafl (2005) 
for a sample of groups (kr < 2 keV). Our sl ightly higher 
slope, compared to that of [Piffareti et alJ J2005t) . can be ex- 
plained by a slight difference in the temperature profile shape. 
We note simply that the samples are not equivalent and that a 
detailed investigation of the slope of the entropy profile will 
require a carefully controlled, unbiased sample spanning the 
mass range from groups to massive clusters. Irrespective of 
the exact slope, that we find S SSr 1 " 1 is interesting considering 
that this va lue is expected for shoc k heating in spherical col- 
lapse (e.g.. lTozzi & Normarl l200lT) . This is further discussed 
in Sec.|5] 
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Fig. 5. Scaled entropy profiles. The radius is expressed as a function 
of density contrast, 6. The entropy has been scaled with the best-fitting 
entropy-mass scaling S2500 x h(z)~ 213 M®-* 4 , using the measured total 
mass, M 2 oo, in units of 10 15 M Q . The shaded grey area corresponds to 
the region enclosed by the mean plus/minus the lcr standard deviation 
(including MKW9). The solid line shows the power law relation, S oc 
(T 66 , obtained from fitting points with 8 < 15000 excluding MKW9. 
The dashed line is an extrapolation of this best fit. 

4.2. Scaled entropy versus 5 

In Figure|5J we show the entropy profiles scaled using the rela- 
tion 5 2500 K h{z} 2 ^ M®-^, where is the total mass M200 in 
units of 10 15 M . Radii are plotted in terms of the density con- 
trast, 6. The shaded area in Fig.[5]shows the region enclosed by 
the mean plus/minus the 1 cr standard deviation. 

Again, while we observe similarity in the entropy profiles 
outside the core regions (6 < 15000), there is a significant in- 
crease in the dispersion at higher density contrast. It is also 
clear that MKW9 is a distinct outlier. This may indicate that 
the mass of this cluster is underestimated, or its temperature is 
overestimated. Excluding MKW9 and fitting the scaled entropy 
profiles with a power law using the BCES method, the slope 
varies from -0.76 + 0.03 when fitting 5 < 3 x 10 5 (~ 0.01 R m ), 
to -0.66 + 0.03, when fitting 5 < 15000 (~ 0.1 R 200 )- Thus the 
entropy does not strictly behave as a power-law in density con- 
trast. The intrinsic dispersion drops from 30 per cent to 10 per 
cent depending on the range of density contrast fitted. Outside 
the core regions, the intrinsic dispersion is smaller than for 
temperature-scaled profiles, in agreement with the lower dis- 
persion observed for the S-M relation compared to that of the 
S—T relation. 



5. Discussion 

Our results put into evidence two main departures from the 
standard self-similar model of cluster formation. First, beyond 



the core region (r > 0.1/? 200) the entropy profiles obey self- 
similarity, having a shape consistent with expectations, but with 
a modified temperature (or mass) scaling. The scaling relations 
are shallower than expected. Second, there is a break of similar- 
ity in the core region: the dispersion in scaled profiles increases 
with decreasing radius. In this section, we will discuss in turn 
both characteristics. 

5.1. Entropy normalisation 

The modified entropy scaling indicates that there is an excess 
of entropy, in low mass objects relative to more massive sys- 
tems, as compared to the expectation from pure shock heating. 
A comparison with adiabatic numerical simulations allows us 
to quantify the absolute value of the excess and to examine 
whether an entropy excess is also present for the most massive 
systems. 

shows the results of adiabatic numerical sim- 
ulations of 30 clusters spanning a mass range of more than a 
factor of ten. Once scaled by the characteristic entropy of the 
halo, 



5200 - \ 
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15 f b H(z) 
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where fi, is the baryon fraction, the simulated profiles are 
closely self-similar, and can be well fitted in the radial range 
~ 0.1-1.0 ^200 by the power-law relation S/S 200 = 
1.26(.R/2?200) U . Assuming f b = 0.14 (a b h 2 = 0.02 and Q m = 
0.3) and typical elemental abundances, we can scale our ob- 
served entropy values to S 200 using the expression 
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2471 keVcm 2 /\l x 15M G 
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Further scaling the radius by the measured value of R^oo allows 
us to compare our results directly with the adiabatic simula- 
tions. In the left hand panel of Figure |6j our observed en tropy 
profiles are compared with the prediction of lVoill J2005D. The 
right hand panel shows the ratio between the best-fitting adia- 
batic power-law relation and the observed data. 

It can be seen that the richer systems all have entropies 
which are in good agreement (both in slope and normalisa- 
tion) with the adiabatic prediction, denoted by the solid line 
in Fig.|6] On average, their entropy is only slightly higher than 
predicted (by ~ 20 per cent), although the effect is not very 
significant. We recall that there is also a ~ 30 per cent differ- 
ence in normalisation between the obser ved M-T relation an d 
that predicted by adiabatic simulations dArnaud et all 12005). 
Interestingly, this corresponds to a ~ 20 per cent entropy ex- 
cess at a given mass for T oc M 2 ^ . The (slight) excess of en- 
tropy in massive systems is thus consistent with a simple in- 
crease of the mean temperature, i.e., of the internal energy of 
the ICM. However, Fig.[6]shows explicitly that the poorer sys- 
tems have a systematically higher entropy normalisation than 
the richer systems. There is approximately 2.5 times more en- 
tropy at 0.2/?2oo in the ICM of A1983, the poorest cluster in 
our sample, than that predicted by gravitational heating. This 
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Fig. 6. A comparison of our observed entropy profiles with the prediction from the adiabatic numerical simulations of TVoitl Left panel: 

The observed entropy profiles have been scaled to Sjoo using Equation 4 The solid line represents the best-fitting power law relation found by 
IVoitl |2005) from fitting adiabatic numerical simulations of 30 clusters in the radial range 0.1 < R 2 oo < 1-0. Right panel: Ratio between the 
observed profiles and the best-fitting power law relation found by Voit 1 2005 ). The dashed line shows the same relation multiplied by 1.2 to 
take into account the difference of 30 percent between the observed normalistion of the M-T relation and that found in adiabatic simulations. 



excess shows that the density of the ICM is also affected at 
lower mass. 

IVoit et all d2003l) and IPonman et all J2003I) independently 
noted that the ICM entropy is highly sensitive to the density 
of the incoming gas and suggested that a smoothing of the gas 
density due to pre-heating in filaments and/or infalling groups 
would boost the entropy production at the accretion shock. The 
entropy amplification is more efficient in low mass systems, 
because they accrete smaller halos more affected by smoothing 
due to pre-heating. In this interpretation, no substantial isen- 
tropic core will develop because the amount of initial preheat- 
ing needed is substantially less than the characteristic entropy 
of the final halo. Clusters will thus be self-similar down to low 
mass but with a modified normalisation scaling as discussed 
above. This is in agreement with our results outside the core 
region (r > 0.1fi 2 oo)- 

Recent numerical simulations which mimic preheating 
by imposing a minimum entropy floor at high z have con- 
firme d the entropy amplifi cation effect due to smooth accre- 
tion (Bor gani et all 120051) . However, the effect seems to be 
substantially reduced when cooling is also taken into account. 
Furthermore, the physical origin of the preheating is still un- 
clear. Heating by SNe seems to be too localise d to have a sig- 
nifica nt effect in smoothing the accreting gas dBorgani et all 
120051) . While AGN might be better candidates to produce 
the extra heating, the observed normalisations require that the 
AGN affects the entropy distribution at least up to riooo- Some 
rece nt theoretical investigations suggest that this may be possi- 
ble jRovchowdhurv et allEooltlVoit & Donahue! 120051) . 



5.2. Localised modification 

While filamentary pre-heating may explain the increased nor- 
malisation of poor systems relative to hotter systems, it does 
not explain the increasing scatter towards the central regions 
in scaled profi les (Figs. El |5] and |6j. The adiabatic numerical 
simulations of Voit (2005) show both a flattening of the slope 
and an increase in the dispersion of the scaled entropy profiles 
in the central regions (< 0.1 Rioo)- However, the dispersion in 
our observed profiles (~ 60 per cent) greatly exceeds that of the 
simulations (~ 30 per cent, cf Fig.[6]and Fig. 11 of lVoil20 05). 
Six clusters out of our total sample of ten (A1991, A2597, 
A 1068, A478, PKS0745 and A2204) have remarkably similar 
scaled entropy profiles, displaying power-law behaviour down 
to the smallest radii measured. Fitting these scaled profiles in 
the radial range 0.01 -0.1 R 200 with power-law using the BCES 
method, we find a slope of 1 . 1 3 +0.05, with a dispersion of only 
1 3 per cent. This slope is very similar to that found by fitting the 
r > 0.1 /?200 range. These six clusters all appear to host a bona 
fide cooling core, each having a centra l temperature decrement 
of a factor ~ 2 (cf lKaastra et all 12004). Similar power-law pro- 
files were found by i n a sample of 13 cooling core clusters by 
Piffa reti et all J2005I) . Strong radiative cooling thus appears to 
generate entropy profiles which display power-law behaviour 
down to very small radii (Fig. [6}. This is simil ar to the quasi- 
steady -state entropy profile in the models of iMcCarthv et al] 
(2004), which include radiative cooling. 

The other four clusters in our sample are characterised by 
a smaller central temperature decrement, larger cooling times 
and shallower entropy profiles. Clearly, some mechanism has 
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modified the cooling history of these clusters. Energy input 
from AGN is regularly invoked as a way of moderating cooling 
at the centres of galaxy clusters (e.g., IChurazov et all l200ll 
and references therein). Our sample contains four clusters 
which have X-ray evidence for interactions between radio and 
X-ray plasma (A478, A2204, A2597, and PKS0745), and yet 
the entropy profiles of all of these clusters increase monotoni- 
cally outward in the canonical fashion. However, we note that 
Chandra observations of clusters with moderate to strong radio 
sources (including A478), show some evidence for a slightly 
shallower slope (~ 0.8 - 0.9 ) but only in the very inner re- 
gions (R < 0.05 R?m\ see e.g..lBe lsole. Birkinshaw & Worrall, 
20051: LFohnstone et all 12005 : San derson. Finoguenov & M ohr. 
2004). The effect of AGN activity on the entropy profile will 
depend on whether the energy input is catastrophic (i.e., oc- 
curring in strong bursts) or distributed (more moderated input). 
It is clear that if the heating is catastrophic in nature, no clus- 
ter has yet been found in which this is evid ent at least from 
the point of view of the en tropy (although see McNamara et all 
2005i iNulsen et all 120051) . Th e heating is likely more dis- 
tributed, via e.g., weak shocks (Fabi an et all r2003). thus pre- 
serving the generally increasing form of the entropy profile. 

Merging events ca n result in substantial mix ing of high and 
low entropy gas (e.g. iRicker & Sarazi rl l200ll) . Relatively lit- 
tle attention has been given to this pos sibility in the literature 
(see, h owever. Fig . 5 oflVoit et ali 2003. Fig. 12 of Belsole et all 
l2004l and Fig. 6 of lPratt. Bohringer & Fino guenovl2005l) . Such 
redistribution of entropy will depend on the scale of the merger, 
whether the merger has disrupted the structure of the cool 
core, and the timescale for re-establishment of the cool core 
if disrupted. In the current sample MKW9, A2717, A1413 and 
A1983 all have flatter core entropy profiles. We note that the 
morphological information for the present samp le would argue 
against recent merger activity in these clusters dPratt & Ar naud 
2005; Papers I and II). However, this does not rule out en- 
tropy modification due to a more ancient merger, particularly 
if the relaxation tim e is less than the cooling time (see also 
lBelsoleetalll2004 . 

6. Conclusions 

The present sample of ten relaxed clusters constitutes the first 
with precise mass data and wide temperature coverage, allow- 
ing the detailed study of the scaling of ICM entropy with both 
temperature and mass. The entropy profiles are sampled with 
good spatial resolution up to ~ 0.5 R200, allowing us to exam- 
ine structural properties avoiding the stacking analysis and ex- 
trapolation on which previous ROSAT/ASCA works relied. 

We have found that the entropy profiles of the present clus- 
ter sample are self-similar beyond ~ O.I.R200, with a shape 
not significantly different from that expected from shock heat- 
ing. The entropy scaling relations are shallower than expected 
from standard self-similar models and adiabatic numerical sim- 
ulations. The slopes of the entropy scaling relations are inde- 
pendent of radius, reflecting the structural self-similarity. The 
entropy-mass (S-M) relation is consistent with the observed 
entropy-temperature (S-T) and mass-temperature (M-T) re- 
lations. The dispersion is smaller about the S-M relation, re- 



flecting the fact that the mass of a cluster is its most funda- 
mental characteris tic. These results co nfirm the trends seen in 
smaller samples (Pratt & Arnaudll2005l). and in stacking anal- 
ysis of ROSAT/ASCA data dPonman et alll2003l) . 

Comparison of our observed profiles with adiabatic numer- 
ical simulations has allowed us to quantify the entropy excess 
relative to that expected in pure gravitational collapse. The ex- 
cess is only ~ 20 per cent in our highest-mass systems. This 
is in line with a simple increase in the mean temperature con- 
sistent with the observed difference in normalisation between 
the observed M-T relation and that predicted by adiabatic sim- 
ulations. However, the excess in low mass systems can be of 
more than a factor of 2.5, indicating that the density of the ICM 
is also affected in these systems. These characteristics are in 
qualitative agreement with a scenario in which entropy produc- 
tion is boosted at the accretion shock. A plausible candidate 
mechanism is smoothing of the accret ed gas due to preheating 
dVoit et allEooilPonman et alll2003f) . 

We have found new and compelling evidence for an in- 
crease in the entropy dispersion in the core regions ~ 60 per 
cent at 0.027?2oo- However the five cooling core clusters in our 
sample have remarkably self-similar power-law profiles with a 
dispersion of only 13 per cent between 0.01 and 0.1 ^200- The 
observed increase in dispersion towards the central regions ar- 
gues for localised entropy modification mechanisms. We con- 
clude that AGN activity and/or transient entropy modification 
due to merging events are good candidates for the modification 
of cooling in these clusters. 

These will be interesting questions to address with numer- 
ical simulations. On the observational side, a larger, unbiased, 
sample of clusters will undoubtedly provide greater insights. 
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Appendix A: Surface brightness profile fit results 

The co-added EPIC azimuthally averaged, vignetting cor- 
rected, background subtracted surface brightness profile was 
computed in the [0.3-3.] keV band for each cluster. The 
profile was corrected for radial variations in the emissivity 
due to abundance or temperature gradients as described in 
Pratt & Arnau d (2003). The gas density profile was obtained 
from fitting parametric analytic models, convolved with the 
EPIC PSF, to the corrected profile. As discussed in Paper I, 
we tried various models and empirically chose the best fitting 
model using the x 2 statistic as a measure of the goodness of fit. 
The models used were : (i) the double /3 (BB) density model 
dPratt & Arnaul l2002l Eq. 3-5); (ii) the modified double /3 
(KBB) model, allowing a more concentrated gas d ensity dis- 
tribution towards the centre (Pratt & Arnaud, 2002, Eq. 6-8), 
and (iii) the sum of three /3 (3B) models for the emission mea- 
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Table 4. Columns: (1) Cluster name; (2) Parametric model; (3) Central density of principal fi model (xl0~ 2 cm" 3 ); (4,5) Central density 
of second and third j3 model components (xlCT 2 h 1 ^ cm" 3 ) (6) Core radius of inner y6-model (arcmin); (7) Core radius of principal /3-model 
(arcmin); (8,9) Core radius of second and third /? model components (arcmin); (10) /? parameter; (11) Radius of boundary between /3 model 
components; (12) £ parameter; (13)^ 2 /d.o.f. 
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sure, in which a comm on value of $ is used to ens ure smooth 
behaviour at large radii ( .Pointecouteau et all l2004i Eq. 2-3). 
The BB and KBB models can be written: 



r < r cut n e {r) = n e] 
r > r cut n e (r) = N 



^(0 



(A.l) 



where £ = 1 for the BB model and £ < 1 for the KBB model, 
and 



N = n e 



with 



An =fi 



2£ 



(A.2) 



The density profile of the BBB model can be written: 



(A.3) 



n e (r) = Jn\ (r) + n\ (r) + n\ (r), 
with 



n ei (r) = n Ci 



1 + 



(A.4) 



(A.5) 



The best-fitting models, and model parameters, are listed in 
Table|H 
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